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Abstract 

With construction halfway complete, IceCube is already the most sensitive neutrino telescope ever built. A rearrange- 
ment of the final holes of IceCube with increased spacing has been discussed recently to optimize the high energy 
sensitivity of the detector. Extending this baseline with radio and acoustic instrumentation in the same holes could 
further improve the high energy response. The goal would be both to detect events and to act as a pathfinder for 
hybrid detection, towards a possible larger hybrid array. Simulation results for such an array are presented here. 
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1. Introduction 

Detection of astrophysical neutrinos at energies 
above 10^^ eV would test cosmic ray production 
models and fundamental particle physics at ^100 
TeV centre of mass energy. 

Techniques beyond optical measurements are re- 
quired to detect 10 or more events per year with a 
reasonable amount of instrumentation. Integrating 
both radio and acoustic receivers with optical instru- 
mentation could provide the necessary sensitivity as 
well as allow for cross calibration and background 
rejection with hybrid events. The full simulation of 
such a detector, combining optical, radio and acous- 
tic detection methods, was done in [l] and revealed 
the possibility to detect ^ 20 GZK neutrinos per 
year, 40% of which were in the form of hybrid events 
detected by more than one technique. 



The optical Cherenkov detection technique is well 
established and fully exploited in IceCube, located 
in the deep ice at the South Pole and currently 
halfway constructed [2 . The radio and acoustic de- 
tection methods are currently in an intensive R&D 
phase, with several experiments (RICE j^J, AURA 
0] and SPATS [5] [6]) underway to measure in situ 
the parameters necessary to establish the optimal 
geometry for a large detector. 

A first step towards the construction of a multi- 
km"^ hybrid detector could consist in a high energy 
hybrid extension around IccCube, combining radio 
and acoustic sensors with optical instrumentation. 
A simulation of a detector with such an extension 
has been performed to study what the improvement 
would be at highest energies. 
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2. Simulation of a possible UHE extension 
to IceCube 

2.1. Geometry 

IceCube will be made of 80 strings arrayed hexag- 
onally with 125 m horizontal spacing. Each string 
is equipped with 60 DOMs (Digital Optical Mod- 
ules), with 17 m vertical spacing at depths between 
1.45 km and 2.45 km. In the geometry simulated, 12 
of the strings are rearranged at 1 km distance from 
the centre, taking into account the presence of the 
skiway on one side of the detector, which does not 
allow for drilling [7j. Half of these strings, plus one 
at the centre of the array, are equipped additionally 
with other technologies: 5 radio antennas are placed 
every 100 m from 200 m to 600 m depth, and 60 
acoustic sensors are located every 15 m from 215 m 
to 1100 m depth. An illustration of the geometry is 
shown in Fig. [T] 
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Fig. 1. Geometry considered for the UHE extension to the 
IceCube detector. Simulated instrumentation, as a function 
of the depth along a single string, is schematically indicated 
on the right. 



2.2. Simulation Methods 

The optical simulation is run with standard Ice- 
Cube software which takes into account parti- 
cle propagation for muons and electrons and the 
Landau-Pomeranchuk-Migdal (LPM) effect. Light is 
propagated using an ice model with dust layers [8] . 
The DOMs are simulated as realistically as possible. 
Both a global and a local trigger condition are set, 
requiring a minimum of 8 modules in the whole ar- 
ray being hit in 5 ^s and at least two adjacent DOMs 
being hit in 1 /is. 



For the radio simulation, the signal from a shower 
is modelled as in [10], including the LPM effect as 
prescribed in [llj . with the assumption of a maxi- 
mum field attenuation length of 1.2 km. Every radio 
receiver is a single half- wave dipole antenna with an 
effective height of 0.27 m at the frequency of 200 
MHz. A Gaussian bandwidth with a a ~60 MHz 
and a sharp high pass filter at 110 MHz are included 
in the receiver model. The trigger condition applied 
requires a signal above 4.5 akT in 4 antennas, where 
CTfeT is the rms noise (including thermal and 250 K 
system temperature) . 

The acoustic simulation uses a modified version of 
the Nishimura-Kamata-Greisen (NKG) parameter- 
ization for hadronic showers, which includes length- 
ening by the LPM effect [12] [13]. Electromagnetic 
showers generated from primary or secondary elec- 
tron neutrinos are expected to be strongly elongated 
by the LPM effect, causing a weaker acoustic sig- 
nal; therefore they are ignored. Given the energy 
deposited at each interaction, the bipolar acoustic 
waveform is simulated in a way similar to that used 
for S AUND ITS' . The signal is calculated taking into 
account a depth-dependent attenuation model |14| , 
which predicts an attenuation length between 7.5 
km and 4.5 km at the instrumented depths. A uni- 
form frequency and angular response of the sensor 
is assumed. The trigger condition requires a signal 
above 9 mPa recorded by at least 3 sensors. The lat- 
est measurements made in situ with SPATS com- 
pared to calibration data in the laboratory show that 
this assumption is reasonable [16j . 

3. Simulation results 

3.1. Neutrino Set 

A set of 10^ neutrinos (including neutrinos and 
antineutrinos of all 3 flavours z/o, t'^j, i't), following 
a E-'^ spectrum from lO^''-^ eV to lO^^-^ eV, has 
been produced using the neutrino generator ANIS 
[17j . This includes propagation through air, rock and 
ice, with CTEQ5 cross sections to calculate neu- 
tral and charged current interaction rates. Neutri- 
nos were generated up to a zenith angle of 120° since 
at these energies upgoing neutrinos are absorbed by 
the Earth. 

The generation volume was a cylinder of 10 km 
radius, 20 km length, resulting in a volume of '^6280 
km"^ . This volume is rotated in the direction of each 
neutrino. All the neutrinos reaching the generation 
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volume are forced to interact and a weight is as- 
signed to take into account the interaction probabil- 
ity. Each interaction can produce a signal detectable 
by any of the techniques which are simulated inde- 
pendently. The identification of the event is unique 
so it is possible to identify events detected by any 
combination of methods. 

The number of events detected by each single 
method and by each combination of methods has 
been used to calculate the effective volume as a func- 
tion of energy. The results are shown in Fig. [2l The 
flatness of the optical curve stems from the contri- 
bution of muons to the signal events. At the energy 
studied here, due to the limited size of the simulated 
interaction volume, the likelihood of a muon to reach 
the detector is independent of the energy. Acoustic 
and radio detection methods are instead sensitive to 
signals produced by hadronic showers with a par- 
ticular threshold and a strength increasing with the 
energy. 

The distribution of hybrid events detected by 
more than one method is shown in Fig.[3]and in Fig. 
3] using the IceCube coordinate system, whose ori- 
gin is located in the centre of the detector, 1948.07 
m below the surface. The horizontal coordinates are 
(x, y) and z is positive upward. 
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Fig. 2. Effective volume for the optical, radio and acoustic 
detection method in the energy range of the simulated neu- 
trino sample. 

The resulting effective volumes have been folded 
with the Engel-Seckel-Stanev GZK flux model [l8] . 
assuming a source evolution parameter ^1^=0.7. In- 
teractions by Vr and 777 arc introduced in this model 
assuming a resulting final mixing 1:1:1 due to flavour 
oscillation. 
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Fig. 3. Distribution of the hybrid events in the x, y coordi- 
nates. 
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Fig. 4. Distribution of the hybrid events in the x, z coordi- 
nates. 



The number of GZK neutrinos detected per year 
by each single method and each combination of 
methods is shown in Table [1] The total number has 
been calculated considering all the events detected 
by any method. 

Reasonable background rejection may require a 
cut on the optical channel multiplicity. With NCh > 
100 the number of optically observable events would 
be reduced by less than a factor of ~2, as shown in 
Fig. El 
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Table 1 

GZK rates per year from the simulated geometry. "Optical" 
refers to IceCube plus the optical channel of the extension. 



Detection option 


GZK events per year 


IceCube 


1.78 


Optical 


3.09 


Radio 


1.31 


Acoustic 


0.16 


Optical + Radio 


0.15 


Optical + Acoustic 
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Radio + Acoustic 
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Optical + Radio + Acoustic 


0.01 


Total 


4.32 
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Fig. 5. Rate of events in the optical channel vs. the number 
of channels hit chosen as trigger threshold. 

4. Conclusion 



The simulated extension increases the number of 
GZK neutrino events seen by IceCube by a factor 
of ^2. The detection of a few additional events by 
the radio method would strengthen the observation. 
The detector extension in the present form would 
be too small to profit from acoustic measurements. 
However in the present effort neither the string nor 
the sensor geometry nor the trigger conditions have 
been optimized for the radio and acoustic detector 
components. 

An optimized medium size extension of IceCube 
could nevertheless build the bridge to a large vol- 
ume, high sensitivity GZK neutrino detector. It 
would also be a valuable prototype to test and 
develop the most suitable technology and infras- 



tructure. 
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